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Parallel  channel  flow  field  with  manifold  openings  is  widely  used  in  Proton  exchange  membrane  fuel 
cells  (PEMFCs)  because  of  its  low-pressure  drop  and  easiness  of  manufacture.  This  research  presents 
a  hydrodynamic  model  to  describe  the  airflow  distribution,  and  the  predicted  pressure  differences  are 
validated  by  experiments.  We  also  investigate  the  influences  of  the  flow  rate,  the  geometry  of  header  and 
the  length  ratio  of  manifold  opening  to  header  region  on  the  airflow  distribution.  Therefore,  the  optimal 
strategy  is  proposed  based  on  an  overall  consideration  of  uniformity  and  configuration  in  the  fuel-cell 
plate  for  application. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Proton  exchange  membrane  fuel  cells  (PEMFCs)  are  considered 
as  the  most  promising  energy  conversion  systems  for  vehicles 
and  station  applications  in  the  future  [1].  However,  many  of  the 
fundamental  processes  are  only  partially  understood,  and  the  con¬ 
figuration  of  the  gas-distribution  plate  is  one  of  the  key  issues  that 
impact  the  commercialization  process.  The  flow  field  should  be 
easier  for  gas  transport  with  less  pressure  drop  and  can  maintain 
uniform  flow  across  the  surface  of  the  electrode. 

Recently,  many  channel  configurations  of  distributor  plates  have 
been  studied,  such  as  parallel  [2-4],  serpentine  [5,6]  and  interdig- 
itated  [7]  channels.  As  for  the  parallel-channel  plate,  the  design 
of  header  configuration  is  quite  important.  Former  researchers 
often  concentrate  on  the  Z-  and  U-type  configurations,  which  have 
been  solved  with  hydrodynamic  models  or  CFD  simulations  [8]. 
In  practice,  the  parallel-channel  plates  with  manifold  openings  [9] 
(see  Fig.  1)  are  often  used  to  avoid  overmuch  pressure  drop  and 
obtain  uniform  flow  distribution.  However,  there  are  few  mod¬ 
els  to  describe  the  flow  distribution  in  such  configuration.  As  a 
matter  of  experience,  the  length  ratio  of  manifold  opening  to 
header  region  can  greatly  influence  the  flow  distribution  in  the 
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plate,  thus,  the  reasonable  value  should  be  discussed  to  assist  the 
designer. 

It  is  the  purpose  of  this  paper  to  provide  a  simple  hydrodynamic 
model  for  the  flow  distribution  in  the  parallel-channel  plate  with 
manifold  openings.  The  experimental  studies  are  also  given  to  vali¬ 
date  the  pressure  difference  distribution  in  such  small  size  of  these 
channels.  In  this  way,  the  key  factors  that  affect  the  flow  distribu¬ 
tion  will  be  found.  Especially,  the  length  ratio  of  manifold  opening 
to  header  region  is  considered.  Therefore,  the  optimal  strategy  for 
both  uniform  distribution  and  compact  construction  is  discussed 
below. 

2.  Model  description 

2.1.  Flow  distribution  in  a  gas -distributor  plate 

The  flow  pattern  in  the  parallel-channel  plate  with  manifold 
openings  is  depicted  in  Fig.  2.  All  the  parallel  channels  are  separated 
by  ribs  and  numbered  in  sequence.  There  are  N  discrete  channels 
in  the  plate;  consequently,  the  header  region  can  be  discretized  on 
a  mesh  of  N  nodes,  and  the  manifold  openings  cover  Nb  of  the  total 
N  channels  in  both  inlet  and  outlet.  The  gas  enters  from  the  inlet 
manifold  opening  and  flows  though  individual  channels  via  the  feed 
header  region,  then  converges  at  the  outlet  distribution  region  and 
leaves  the  plate  from  the  outlet  manifold  opening.  The  volume  flow 
rate  in  each  of  N  individual  channels  is  denoted  by  qx  ( i ).  Thus,  in 
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Nomenclature 

List  of  symbols 

A 

cross-section  area  (m2) 

d 

width  (m) 

D 

hydraulic  diameter  (m) 

err 

error 

/ 

friction  factor 

Ry 

force  exerted  on  control  volume  (N) 

h 

height  (m) 

L 

length  (m) 

M 

mass  flow  rate  (kg  s-1 ) 

N 

number  of  channels 

Nb 

number  of  channels  covered  by  the  manifold  open- 

ing 

P 

perimeter  (m) 

Pin 

pressure  at  the  beginning  of  the  channel  (Pa) 

Pout 

pressure  at  the  end  of  the  channel  (Pa) 

Qx 

volume  flow  rate  through  each  channel  (m3  s-1 ) 

a 

cumulative  volume  flow  rate  (m3  s-1 ) 

Ux 

mean  velocity  in  the  channel  (m  s-1 ) 

Greek  letters 

a 

aspect  ratio 

air  viscosity  (Pas) 

P 

air  density  (kgm-3) 

X 

shear  stress  (Nm-2) 

0 

characteristic  parameter  (kg  m-3  s-1 ) 

V 

non-dimensional  geometry  parameter 

Subscripts 

in 

near  the  inlet  side 

out 

near  the  outlet  side 

X 

direction  along  the  channel 

y 

direction  along  the  header 

0 

standard  situations 

the  inlet  distribution  region,  the  cumulative  volume  flow  rate  (Qjn 
(i))  is  defined  as 

i- 1 

Qin(i)  =  J2qxU)  =  Qin(i  -  1 )  +  <?*('  -  1 ).  (i  =  2,  3,  . . . ,  N  -  Nb) 

J  =  1 

(1) 

Qin(l)  =  0  (2) 

In  the  inlet  manifold  opening,  the  mass  balance  equation  can  be 
written  as 

N  N 

Qm  =  Y.qx(j)=  X!  ‘fcW  +  QinUV-Nb)  (3) 

j=l  J=N-Nb 

where  Qin  is  the  total  feed  volume  flow  rate,  while  in  the  outlet 
distribution  region,  the  similar  variables  are  given  by 

N 

Qout(Q  —  ^  ]qx(j)  —  Qout(i  +  1)  +  Qx(i )? 

j=i 


(i  =  Nb  +  l,Nb  +  2,...,N-l) 


(4) 


Fig.  1.  Schematic  diagram  of  parallel-channel  plate  with  manifold  openings  in  the 
header  region.  (1)  Header  region,  (2)  manifold  opening,  (3)  distribution  region,  (4) 
flow  field  region,  (5)  channel  and  (6)  rib. 


In  the  outlet  manifold  opening,  the  mass  balance  equation 
becomes: 

Qout  =  Qx(j )  +  Qout(Nb  +  1 )  (6) 

i=  i 

Because  we  suppose  there  is  neither  electrochemical  reaction  nor 
leakage  in  the  flow  field,  the  Eqs.  (3)  and  (6)  are  equal. 

2.2.  Steady  channel  flow  with  lower  Re 


Now,  the  flow  distribution  in  the  plate  can  be  obtained,  and  the 
flow  rate  through  each  channel  is  governed  by  the  pressure  dif¬ 
ference  between  the  beginning  and  end  of  the  channel,  which  is 
expressed  as  Pin(i)  minus  Pout(i).  For  steady  Hagen-Poiseuille  flow 
in  a  channel,  the  net  pressure  balances  the  force  caused  by  the  shear 
stress  at  the  wall: 


(Pin(Q  —  Pout(0)Ax  —  TxPxLx 


(7) 


where  Ax  is  the  channel  cross-section  area,  Px  is  the  channel  perime¬ 
ter,  and  Lx  is  the  length  of  the  channel.  The  shear  stress  ( rx )  is  usually 
expressed  in  terms  of  the  friction  factor  fx. 


fx  — 


(8) 


1/2  pux(i)2 

where  ux{i)  is  the  mean  velocity  in  the  channel,  and  the  friction 
factor  is  a  function  of  the  Reynolds  number: 

pDxux(i ) 


Rex  = 


P 


(9) 


where  Dx  represents  the  hydraulic  diameter,  and  is  defined  as 


(10) 


By  the  empirical  correlation  of  Kays  and  Crawford  [10],  the  ( Ref)x 
depends  on  the  channel  aspect  ratio  (ax  =  dx/hx),  hence: 


( Ref)x 


13.74+  10.38  exp 


(ID 


From  Eqs.  (7)-(10),  the  pressure  difference  from  the  inlet  to  the 
outlet  is  described  as 


(5)  PintO  —  Pout(0  —  <t’xLxux(i) 


(12) 


Qout(N)  =  <JX(N) 
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Qout  (N) 
Qout  (N-1) 
Qout  (N-2) 

Qout  (N-Nb+1) 
Qout  (N-Nb) 
Qout  (N-Nb-1) 


Qout  (Nb+1) 


Fig.  2.  Illustration  of  the  flow  distribution  and  mesh  network  in  a  fuel-cell  plate  with  manifold  openings. 


here 


<PX  = 


2  (Ref)x 
D2x 


(13) 


Then  the  pressure  difference  from  each  channel  can  be  obtained 
directly  from  the  linear  flow  velocity  (ux(z))  in  the  corresponding 
channel. 


2.3.  Pressure  distribution  in  header  regions 


Derivation  of  the  momentum  equation  is  formulated  for  a  suit¬ 
able  control  volume  located  in  the  outlet  distribution  region,  as 
shown  in  Fig.  3.  For  one-dimensional  axial  flow,  the  momentum 
balance  is  stated  as 


<KpUyUy dv  = 


fy 


(14) 


where  Fy  is  a  force  exerted  on  this  volume  and  is  composed  of  two 
forces:  one  is  due  to  normal  stress  on  the  cross-sectional  flow  area, 
as  the  first  integral  on  the  right-hand  side  in  Eq.  (15);  the  other  is 
derived  from  viscous  drag  at  the  channel  walls.  Restricting  attention 
to  the  steady  flow: 


With  the  help  of  Gauss  divergence  theorem,  the  first  integral  on 
the  right-hand  side  converts  to  a  volume  integral.  The  shear  drag 
term  can  be  rewritten  in  one  dimension: 

[  9(yy)dV=  f  (V  ■  Zy)dAy  -  [  TyPydy  (16) 

Jcv  °y  JCV  Jy 

where  Py  is  the  perimeter  of  the  distribution  region,  and  Ay  is  the 
crossing-section  area.  ry  is  the  normal  stress,  which  depends  on  the 
pressure  drop.  Shrinking  the  control  volume,  the  desired  differen¬ 
tial  equation  can  be  expressed  as 

(J(pUyUy)  dp  Py^ 

(17) 

Similar  to  the  channel  flow,  the  wall  stress  is  determined  in  terms 
of  a  friction  factor  /y,  which  depends  on  the  header  dimensions 
and  fluid  properties.  Assuming  the  properties  are  constant,  for  an 
incompressible  flow,  the  momentum  equation  can  be  rewritten  as 

2pUy<ly  +  +  0yUy  =  0  (18) 


here 


[  d2aldv 

Jcv 


L 


xy  ri  dAy  — 


(15) 


<Py  =  fL 


2  (Ref)y 

D2 

uy 


(19) 


The  subscript  y  =  out  refers  to  the  exhaust  header  region  and  y  =  in 
refers  to  the  feed  header  region. 

Because  there  are  a  discrete  number  of  channels,  it  is  appropriate 
to  represent  and  solve  the  momentum  equations  in  discrete  form. 
As  shown  in  Fig.  3,  the  control  volume  in  the  outlet  distribution 
region  is  taken  as  rectangular,  with  the  height  hy  and  the  width  dy. 
Since  the  differential  equation  is  first-order  and  the  boundary  con¬ 
ditions  should  be  satisfied,  it  is  appropriate  to  forward  difference 
this  equation. 


Pout(i  +  1) —  ^out(0  —  —  2  puout{i){tiout{i  +  1)  —  Uout{i))Sy 


(  2(Ref)y  \ 

Dj  ) 


Uout(i)<5y, 


(f  =  Nb,Nb  +  l,...,N-l) 


(20) 


Fig.  3.  Control  volume  for  a  section  in  the  outlet  distribution  region. 
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where  u0ut(Q  can  be  obtained  from  the  flow  distribution  in  Section 
2.1: 


Uout(Q  — 


Qout(i) 


(1  =  1,  2, 


,N) 


(21) 


Specifically,  for  the  outlet  manifold  opening: 

Pout(Q  =  Pexit,  (*  =  1,  2 . ATb)  (22) 

Similarly,  for  the  inlet  manifold  opening: 

Pin(z)  =  PIN,  (i  =  N-Nb,  N-Nb  +  1,...,N)  (23) 

However,  the  sense  of  the  difference  in  the  feed  header  region  is 
opposite  that  of  the  exhaust  header: 

Pin(Q  —  Pin(*  —  1)  —  — 2pu  m(i  —  \){U[n(i)  —  U[n(i  —  l))<5y 

.2  (Pe/)/ 


-  M 


D2 

uy 


here 


uin(Q  — 


Qin(0 


uin(z-l)5y,  (z  =  2,  3, . . . ,  N  —  Nb)  (24) 


(f=l,2,...,N)  (25) 


2.4.  Solution  algorithm 


The  program  layout  is  shown  in  Fig.  4,  and  this  calculation  starts 
with  an  initially  guessed  set  of  channel  volume  flow  rates,  qx(i). 
Then  the  cumulative  flow  rates  in  the  headers,  Qm(z)  and  Qout(Q, 
iiout(z)  and  uin(i),  can  be  obtained  from  the  Eqs.  (1),  (4),  (21),  and 
(25).  These  velocities  are  used  to  calculate  the  pressures  in  header 
regions.  The  next  step  is  to  calculate  the  corresponding  flow  rates 
in  the  channel,  ux_new(0  and  q*_new(i),  from  the  Eqs.  (12)  and  (13). 
Of  course,  the  new  value  of  flow  rates  should  be  adjusted  first. 

4*dnew(i)  —  N  - -Qx_new{i)  (26) 

/  w_i  Qx_new(0 


At  last,  the  accumulate  error  has  to  be  compared  with  the  set  pre¬ 
cision. 


E" 

2-^i=  1 


Aw©  -  q*(0 

Qx(i ) 


(27) 


The  program  will  keep  on  running  with  a  new  set  of  channel  flow 
rates  until  all  the  flow  rates  are  close  to  the  old  ones  within  the 
error  limit. 


Fig.  4.  The  program  layout  for  flow  distribution  calculation. 


3.  Experimental 

In  order  to  validate  the  flow  distribution  in  the  parallel-channel 
plate,  a  fine  facility  was  used  as  shown  in  Fig.  5.  The  transparent  sep¬ 
arator  was  made  of  PMMA  (polymethyl  methacrylate),  of  which  62 
pairs  of  little  holes  were  punched  out  along  the  edge  in  two  rows. 
The  controlled  airflow  entered  from  the  inlet  manifold  opening, 
and  this  certain  amount  flow  was  achieved  by  the  mass  flow  con¬ 
troller  (D07-19A).  Also  the  pressure  drop  along  each  channel  could 
be  measured  from  the  corresponding  pair  of  holes  in  the  separa¬ 
tor.  In  order  to  avoid  the  edge  effect,  these  holes  were  kept  a  certain 
distance  away  from  the  inlet  or  outlet.  So  the  test  length  of  the  chan¬ 
nel  was  275  mm  of  the  whole  358  mm,  and  the  experimental  data 
should  be  enlarged  proportionally.  Two  pressure  transducers  were 
used  for  pressure  measurements  over  the  range  of  0-3  kPa  (STD110 
Honeywell)  and  0-20  kPa  (STD930  Honeywell)  with  0.065%  full- 
scale  accuracy  for  both.  In  this  way,  any  pressure  difference  between 
the  beginning  and  end  of  the  channel  would  be  received.  The  oper¬ 
ation  was  kept  in  room  temperature  and  other  parameters  of  the 
set-up  were  described  in  Table  1. 


© 


Fig.  5.  Illustration  of  the  facility  with  62  pairs  of  holes  for  pressure  difference  mea¬ 
surement  in  each  channel. 
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Table  1 

Specification  of  the  flow  field  configuration  and  operating  conditions. 


Variable 

Units 

Value 

Number  of  channels 

- 

62 

Total  length  of  the  header  region 

mm 

90 

Length  of  the  channel 

mm 

356 

Room  temperature 

°C 

20 

Air  density 

kg  ITT  3 

1.293 

Air  viscosity 

Pas 

1.72  x  10“5 

Standard  mass  flow  rate 

kgs-1 

2.16  x  10“5 

Table  2 

Experimental  cases  and  the  corresponding  parameters. 

Case  description 

Volume  flow  rate 

Total  pressure 

Vx 

'Ey 

(seem) 

drop  (kPa) 

#1 

1000 

5.77 

8.68 

5.51 

#2 

4000 

21.0 

8.68 

5.51 

#3 

6000 

31.0 

8.68 

5.51 

#4 

1000 

2.59 

7.97 

3.78 

4.  Results  and  discussion 


ro 

CL 

0)  ^ 
3  8 

</>  c 
0  0 

□I 


^>o  -  o  o  o  o  6  o  Q  Q  .  .Q  o  .  .  .0  Q 

Case  1  predicted 
□  Casel  measured 
Case2  predicted 
Case2  measured 
Case3  predicted 
$  Case3  measured 
Case4  predicted 
O  Case4  measured 

■  (a) 

4  4-4-  +  +  4  ■*•  •+  4  -f  4  4  4  4 

D  O  •  Q  Q  Q  □  0  .[}•□ 

0  10  20  30  40  50  60 

Channel  Number 


Fig.  7.  (a)  Compare  the  effect  of  different  flow  rates  on  the  pressure  difference  dis¬ 
tribution  by  model  prediction  and  experimental  measurement,  and  then  (b)  show 
the  corresponding  ratio  of  flow  rate  distribution  map. 


4.1.  Pressure  differences  and  flow  distribution 

In  this  study,  four  experiments  under  different  conditions  are 
listed  in  Table  2.  In  case  4,  the  volume  flow  rate  of  the  air  was 
1000  seem,  with  the  inlet  and  outlet  pressure  2.83  and  0.24  kPa, 
respectively.  Comparing  with  the  model  prediction  in  Fig.  6(a),  the 
measured  pressure  differences  in  the  parallel  plate  are  indicated  as 
stars.  From  Eq.  (12),  the  flow  rate  has  a  liner  relationship  with  the 
pressure  difference  in  each  channel.  Thus,  the  corresponding  flow 
distribution  map  is  drawn  in  Fig.  6(b).  Along  the  numbered  chan¬ 
nel,  the  maximum  volume  flow  rate  is  acquired  at  the  29th  channel, 
which  is  16.6  seem.  By  contrast,  the  minimum  value  lies  at  the  edge 
of  the  gas-distributor  plate,  which  is  15.9  seem  in  the  1st  and  Nth 
channel. 

4.2.  Effect  of  feed  flow  rate  on  flow  distribution 


of  pressure  difference  are  plotted  against  the  number  of  channels 
for  various  feed  flow  rates,  and  they  were  also  validated  by  experi¬ 
ments.  It  is  found  that  the  non-uniformity  in  the  flow  distribution 
increases  with  the  flow  rate.  In  this  configuration,  the  minimum 
flow  rate  lies  at  the  edge  of  the  plate,  and  the  highest  mass  flow 
usually  exists  near  the  centre  (see  Fig.  7(b)). 


4.3.  Two  characteristic  parameters 


With  the  objective  of  generalizing  the  results,  the  characteristic 
parameters  in  Eqs.  (13)  and  (19)  can  be  recast  in  terms  of  non- 
dimensional  parameters: 


&X  = 


&XLXAX 

M0 


%  = 


OyLyAy 

Mo 


(28) 

(29) 


For  a  fixed  fuel  cell  (Cases  1-3  in  Table  2),  the  feed  volume  flow 
rate  determined  the  total  inlet  and  outlet  pressure,  which  could  be 
measured  in  situ.  In  this  experimental  set-up,  the  stable  pressure 
differences  of  all  the  channels  could  be  detected  and  compared  with 
the  model’s  prediction.  As  shown  in  Fig.  7(a),  the  predicted  curves 


.9  S. 

E  £ 
c  0 


„  </) 
</>  0 
CD  0 

o  fr 


■  . 

- Pin  predicted 

- Pout  predicted 

A  P  predicted 
$  A  P  measured 

(a) 

_ 1 _ 1 _ 1 _ 1 _ 1 - 1 - 

0  10  20  30  40  50  60 

Channel  Number 


Fig.  6.  (a)  Validate  the  pressure  difference  distribution  and  (b)  show  the  correspond¬ 
ing  flow  rate  distribution  in  the  parallel-channel  plate. 


As  for  the  application,  it  is  appropriate  to  specify  the  standard 
mass-flow  rate  M0.  Thus  the  non-dimensional  group  ^ *  depends 
on  the  channel  geometry  and  depends  on  the  header  geom¬ 
etry,  without  considering  the  flow  rates.  In  cases  1  and  4,  the 
feed  flow  rates  were  same,  while  the  dimensions  of  the  header 
and  channel  were  different,  which  would  have  an  effect  on  the 
two  non-dimensional  parameters  and  pressure  drop  (see  Table  2). 
By  increasing  the  depth  of  the  channel  and  depth  of  header 
from  0.35  to  0.55  mm  in  our  experiments,  the  total  pressure  drop 
decreased  from  5.77  to  2.59  kPa,  the  channel  geometry  parameter 
(tf'x)  decreased  from  8.68  to  7.97,  and  the  header  geometry  parame¬ 
ter  (Wy )  decreased  from  5.51  to  3.78.  In  order  to  compare  the  effect 
of  the  two  parameters  on  flow  distribution,  the  calculation  results 
are  shown  in  Fig.  8.  It  is  apparent  that  the  header  geometry  param¬ 
eter  (^y)  greatly  affects  the  flow  distribution  in  the  plate,  which 
will  become  more  uniform  when  & y  decreases.  By  contrast,  the 
channel  geometry  parameter  (<^x)  nearly  has  no  effect  on  the  flow 
distribution. 

4.4.  Effect  of  the  length  ratio  of  manifold  opening  to  header  region 

In  our  experiments,  the  manifold  opening  length  (Lb)  is  40  cm, 
while  the  total  length  of  the  header  region  (Ly)  is  90  cm,  thus,  the 
length  ratio  of  manifold  opening  to  header  region  ( a  =  Lb/Ly )  is  4/9. 
To  explore  this  ratio  effect  on  flow  distribution,  different  designs 
are  shown  in  Fig.  9.  When  a  is  increased  from  4/9  (see  Fig.  9(a)),  the 
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Fig.  8.  Compare  the  effect  of  two  geometrical  parameters  on  the  flow  distribution 
respectively,  when  (a)  4fx  =  8.68  and  (b)  4fy  =  5.51. 


Fig.  10.  Comparison  of  strategy  for  uniform  flow  distribution:  (a)  the  geometry 
aspect  effect  and  (b)  the  suitable  length  ratio  of  manifold  opening  to  header  region. 


a  =  8/9 


a  =  6/9 


a  =  4/9 


Fig.  9.  Effect  of  the  length  ratio  of  manifold  opening  to  header  region  on  flow  dis¬ 
tribution:  (a)  a  >  4/9  and  (b)  a  >  4/9. 


opening  to  header  region  is  also  important  for  flow  distribution. 
And  the  opening  length,  which  takes  one  third  of  header  region, 
will  be  fit  for  both  uniform  distribution  and  compact  construction 
in  the  fuel-cell  plate  for  application. 

5.  Conclusions 

The  flow  distribution  in  a  parallel-channel  plate  with  the  man¬ 
ifold  openings  is  analyzed  by  the  hydrodynamic  model,  which  is 
also  validated  by  experiments.  It’s  found  that  the  non-uniformity 
in  the  flow  distribution  increases  with  the  feed  flow  rate.  In  our 
experiments,  the  maximum  flow  rate  usually  exists  near  the  cen¬ 
tre.  The  results  show  that,  the  header  geometry  parameter  (&y) 
usually  has  great  effect  on  the  flow  distribution.  Thus,  in  order  to 
get  uniform  distribution,  the  large  size  and  the  square  cross  section 
of  the  distribution  region  should  be  commended.  In  addition,  the 
effect  of  the  length  ratio  of  manifold  opening  to  header  region  is 
also  discussed,  and  1:3  is  considered  as  the  suitable  ratio  for  both 
uniform  distribution  and  compact  construction  in  the  fuel-cell  plate 
for  application. 


channel-flow  distribution  will  become  more  uniform  and  the  peak 
value  of  flow  flux  usually  exists  in  the  center.  What  is  more  inter¬ 
esting,  the  peak  will  break  up  to  two  sharp  peaks,  when  the  ratio  is 
decreased  from  4/9  (see  Fig.  9(b)).  In  addition,  the  two  peaks  will 
move  from  the  centre  to  the  edge,  respectively,  and  the  center  will 
have  the  lowest  volume  flow  rate.  Apparently,  the  flow  distribution 
will  become  more  uniform  if  this  ratio  is  higher,  however,  it’s  not  fit 
for  practical  application.  Thus  there  is  a  contradiction  for  the  flow 
field  design,  which  needs  to  be  optimized. 

4.5.  Optimized  design  for  uniform  distribution 

By  analyzing  the  Eqs.  (19)  and  (29),  the  hydraulic  diameter  (Dy) 
also  has  effect  on  the  header  geometry  parameter  (^y),  and  it  will 
influence  the  flow  distribution  greatly  as  shown  in  Fig.  10(a).  Simi¬ 
lar  to  Eq.  (10),  the  largest  Dy  will  be  obtained  in  a  square  header,  if 
the  perimeter  is  fixed.  In  general,  the  large  size  and  the  square  cross 
section  of  the  distribution  region  will  be  better  for  uniform  distribu¬ 
tion.  As  illustrated  in  Fig.  10(b),  the  proper  length  ratio  of  manifold 
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